Mon. Not. R. Astron. Soc. 000, [T]-?? (2002) Printed 16 March 2010 (MN WT^X style file v2.2) 



Proper motion and apparent contraction in J0650+6001 



M. D. Dallacasa^'^ 

^ Dipartimento di Astronomia, Universitd di Bologna, via Ranzani 1, 1-40127, Bologna, Italy 

^ Istituto di Radioastronomia - INAF, Via P. Gobetti 101, 1-40129 Bologna, Italy 

^ Instituto de Astrofisica de Canarias, c/ Via Ldctea s/n, E-S8205 La Laguna (Tenerife), Spain 



Received 16 March 2010; accepted ? 



ABSTRACT 

We present a multi-epoch and multi-frequency VLBI study of the compact radio source 
J0650-f 6001. In VLBI images the source is resolved into three components. The central 
component shows a flat spectrum, suggesting the presence of the core, while the two 
outer regions, with a steeper spectral index, display a highly asymmetric flux density. 
The time baseline of the observations considered to derive the source expansion covers 
about 15 years. During this time interval, the distance between the two outer com- 
ponents has increased by 0.28±0.13 mas, that corresponds to an apparent separation 
velocity of 0.39c ± 0.18c and a kinematic age of 360±170 years. On the other hand, a 
multi-epoch monitoring of the separation between the central and the southern com- 
ponents points out an apparent contraction of about 0.29 ±0.02 mas, corresponding to 
an apparent contraction velocity of 0.37c ± 0.02c. Assuming that the radio structure 
is intrinsically symmetric, the high flux density ratio between the outer components 
can be explained in terms of Doppler beaming effects where the mildly relativistic jets 
are separating with an intrinsic velocity of 0.43c ± 0.04c at an angle between 12° and 
28° to the line of sight. In this context, the apparent contraction may be interpreted 
as a knot in the jet that is moving towards the southern component with an intrinsic 
velocity of 0.66c ± 0.03c, and its flux density is boosted by a Doppler factor of 2.0. 

Key words: radio continuum: general - radiation mechanisms: non-thermal - quasars: 
individual: J0650-F6001 



1 INTRODUCTION 

' The evolutionary stage of a powerful radio source originated 
by an active galactic nucleus (AGN) is related to its linear 
size. In this context, compact symmetric objects (CSO), 
which are powerful (Li.4ghz > 10^^ W/Hz) and intrinsi- 
cally small (i.e. linear size LS < 1 kpc) radio sources, should 
represent a young stage in the individual radio source life. 
These objects usually have a roughly symmetric structure, 
with core, jets, and mini-lobes/hotspots, resembling a 
scaled-down version of the large classical radio galaxies 
which conversely have linear sizes of hundreds of kpc 
or even a few Mpc. The main characteristic displayed 
by young radio sources is the convex synchrotron radio 
spectrum that peaks at frequencies in the GHz regime. 
The spectral turn over is usually due to synchrotron self- 
absor ption (SSA; lOrienti fc Dallacasal l2008l : I Snellen et al.l 
I2OOO !). although an additional contribution from free-free 
absorption ( FFA) has been found in the most compact 
objects (e.g. iDallacasa fc Orientil l2009l : [Mutoh et al. I l2002l : 
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iKameno et"ai]|2000l '). 

Conclusive evidence of the genuine youth of this class of 
objects came from the determination of both kinematic and 
radiative ages. 

The determination of the kinematic age is based on 
multi-epoch milliarcsecond-scale resolution observations 
spanning several years, and able to measure at which rate 
the hotspots are increasing their separation. Assuming that 
the separation velocity «sep has maintained constant, the 
kinematic age tkin can be estimated: 

^kin ~ LS ■ t'scp • 

From the an alysis of a dozen of the m ost compact CSOs 
(LS < 20 pc. lPolatidis fc Conwavl[20ol ') . it has been derived 
that the separation speed is generally in the range of 0.1c 
and 0.4c, leading to kinematic ages of a few thousand years. 
The radiative age can be estimated by multi-frequency 
observations able to constrain the radio spectral shape and 
thus to determine at which frequency the spectral break 
occurs. The break frequency z^br is strictly related to the 
radiative lifetime of the synchrotron emitting electrons tsy-n ■ 
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Table 1. Log of the archival VLB A observations analyzed in this 
paper. 



Freq. Obs. date Beam Mfmax Antennas 

GHz mas2 MA 



5.0 22 Nov 1999 2.46x1.94 125 VLBA -Sc -NI 
5.0 17 Dec 2004 1.65x0.88 172 VLBA + Eff 
8.4 17 Dec 2004 0.90x0.54 290 VLBA + Eff 



Once the magnetic field is known, for example assuming 
minimum total energy content corresponding to equipar- 
tition between p articles and magnetic field {H) energies 
l|Pacliolczvklll970l l. and following some assumptions (jMurgial 
|2003| ). the radiative age can be easily estimated by means of: 

tsyn = const ■ l'l^{^H~^^'^ . 

Several studies aimed at estimating the radiative age in 
compact radio objec t s indicate age s between 10'^ - lO"* years 
l|Murgia et al.l 1 19991 : iMurgial |2003| ). in excellent agreement 
with the kinematic ages. 

Given all the assumptions mentioned earlier, these methods 
do not provide the accurate source age, that can be 
improved either by increasing the time interval spanned 
by the observations or through a better sampling of the 
frequency coverage used to derive the radio spectrum. 
Among all the CSOs studied in this framework, the radio 
source J0650+6001 represents a peculiar case. This source 
is identified with a quasar at redshift z = 0.455, and 
its optical spectrum is characterize d by very weak broa d 
lines and prominent narrow lines (|Stickel fc Kuhilll993h . 
It has a non-aligned triple radio structure with a total 
angular size of 7 mas (~40 pc) , and its total radio spectrum 
turns over around 5.5 GHz jOrienti et al.l [ioOTi '). From a 
previo us study no evidence of proper motion has been 
found (jOrienti fc Dallacasal r2008l \ although the resolution 
at 5 GHz (i.e. the frequency at which several observations 
spanning almost a decade have been performed) was not 
adequate to reliably estimate small changes in the compo- 
nent positions. However, another analysis of the position of 
the source compone nts based on global VLBI observations 
l|Akuior et al.l Il996l ). indicated an apparent contraction 
between the central and the southern components. 
In this paper we report on the results on multi-epoch VLBI 
and 5-GHz VLA data of the compact symmetric object 
J0650-I-6001, and we present an interpretation to explain 
the radio properties of this source. 

Throughout this paper, we assume the following 
cosmology: .ffo = 71 km s~^ Mpc~^, JIm = 0.27, and 
Q.A = 0.73, in a flat Universe. At the redshift of the target 
1 = 5.773 kpc. The spectral index is deflned as S{v) oc i^"" . 
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Figure 1. VLBA images at 5 GHz (top and centre panels), and 
at 8.4 GHz (bottom panel) of the target source J0650-I-6001. 



To determine a possible source growth, wc complemented 
the information available in the literature for the target 
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Table 2. Multi-epoch VLA flux density at 5 GHz of the source 
J0650+6001. Column 1: flux de nsity in mjy; Col. 2: observ - 
ing date; Col. 3 : Ref erence: 1: [uivestad fc JohnstonI lll98ll) : 
2: IPerlevI lll982l}: 3: lO'Dea et al.l dlQQOl'): 4: this work: 5: 



iDallacas a et al.r (|2000|'). 6: iTinti et al.l 1I2OO5I '): 7: lOrienti et al] 



l|2007j); 8: lOrienti fc Dallacasal l|2008l ') 



Flux 



Obs. date 



Rcf. 



906 


2 Fob 1979 


1 


830 


18 Nov 1980 


2 


916 


30 Dec 1984 


3 


940 


21 Jul 1991 


4 


969 


23 Feb 1993 


4 


1058 


24 Sep 1995 


4 


1236 


14 Nov 1998 


5 


1136 


3 Jul 2002 


6 


1150 


14 Sep 2003 


7 


1106 


28 Jan 2004 


7 


1180 


19 Nov 2006 


8 



source J0650-(-6001 with archival VLBA data obtained 
in November 1999 at 5.0 GHz, and VLBA-|-Effelsberg 
observations at 5.0 and 8.4 GHz carried out in December 
2004 (see Table [1]). The data reduction was carried out by 
means of the NRAO AIPS package. The accuracy of the 
amplitude calibration has been estimated to be within 5%. 
The final images were produced after a number of phase 
self-calibration iterations. Amplitude self-calibration was 
applied only once at the end of the process using a solution 
interval of 30 min, i.e. longer than the scan length, to 
remove residual systematic errors and to fine tune the flux 
density scale, but not to force the individual data points to 
follow the model. 

To construct a well-sampled light-curve at 5 GHz, we 
complemented the information already available in the 
literature with archival VLA observations (see Table [2]). We 
retrieved and analyzed only datasets in which the source 
had been observed for at least 5 minutes in the C band. 
Since the peak frequency is close to 5 GHz, observations 
carried out at frequencies in the range of about 4.9 and 5.0 
GHz do not cause significant scatter on the determination 
of the fight curve. Data were reduced following the standard 
procedure implemented in the NRAO AIPS package. 
Uncertainties in the determination of the fiux density are 
dominated by amplitude calibration errors which have 
been found to be within 3% in all the datasets considered 
here. VLA images were obtained after several phase-only 
self-calibration iterations. 



3 RESULT 

3.1 Radio images 

Full resolution VLBI images are presented in Fig. [1] At 
8.4 GHz, in addition to the full resolution image, we also 
produced a low resolution image using the same uv-range, 
image sampling and restoring beam of the image at 5.0 
GHz, in order to produce a spectral index map of the 
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Figure 2. Spectral index distribution across the target source 
J0650-I-6001 superimposed on the 8.4 GHz contours convolved to 
the 5 GHz beam. 



source (Fig. [2]). No VLA image is presented, since the 
source is unresolved on the angular scales typical of this 
interferometer. 

For each VLBI image we provide the observing frequency, 
the restoring beam plotted on the bottom left corner, 
the peak flux density in mjy/beam, and the first contour 
intensity (/. c.) in mJy/beam, that corresponds to 3 times 
the off-source noise level. Contour levels increase by a factor 
2. 

Flux density and angular size of each source sub-component 
have been measured by means of the task JMFIT, which 
performs a 2D-Gaussian fit to the source components on 
the image plane. TVSTAT, which performs an aperture 
integration on a selected region on the image plane, has 
been used to measure the total flux density of the source, 
which results in agreement with the sum of the flux densities 
of the various sub-components. The parameters derived 
from the fitting procedure are reported in Table O 



3.2 Morphology 

When imaged with the high spatial resolution provided 
by VLBA-|-Effelsberg observations, the radio source 
J0650-I-6001 shows a non-aligned triple structure, where 
the two main components, labelled C and S in Fig. 
[T] represent almost the 99% of the total source ra- 
dio emission, i n agreement with previous observations 
lIXu et al.lll995l:lAkuior et al.lll996l : IStangheUini et al.lll999l : 



lOrienti fc Dallacasal l2008t ). Gomponent C, accounting for 
63 and 72 mJy at 5 and 8.4 GHz respectively, possesses a 
flat spectrum (afio ~ —0.05 ± 0.15), suggesting therefore 
the presence of the source core. At 8.4 GHz it is clearly 
resolved in the North direction, indicating the presence of 
a jet-like component (labelled J in Fig. [!};) that accounts 
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Table 3. Angular separation of components N and C with respect 
to S measured on images at 5 GHz. Column 1: angular separation 
between components C and S; Col. 2: position angle of compo- 
nent C with respect to component S; Col. 3: angular separation 
between component N and S; Col. 4: position angle of compo- 
nent N with res pect to compone nt S; Col. 5: observing d ate; Col. 
6: refere nce: 1: IXu et al.l jiggsll; 2: iBritzen et"all |[2to3); 3: this 
work; 4: l0rienti fc Dallacasal l|200Sh . 



dc 



<t>C 
deg 



0N 
deg 



Observing date Ref. 



3.03 
3.08 
2.96 
2.90 
2.81 
2.79 



35.2 
36.4 
36.1 
35.9 
35.3 
35.6 



6.29 23.4 



6.56 20.6 



Jun 1991 
17 Sep 1994 
9 Feb 1998 
22 Nov 1999 
17 Dec 2004 
28 Jul 2006 




year 

Figure 3. The contraction between the central and the southern 
components of the source J0650-I-6001. The line represents the 
best fit of = —0.02 mas/year. 



Table 4. Observational parameters of the components of 
J0650-I-6001. Column 1: source component; Columns 2,3: flux 
density at 5 GHz from observations carried out in 1999 and 2004 
respectively; Col. 4: flux density at 8.4 GHz; Col. 5: spectral in- 
dex between 5 and 8.4 GHz; Cols. 6, 7: deconvolved major and 
minor axis, respectively; Col. 8: position angle of the major axis. 
Spectral index and component angular size are from observations 
carried out in 2004. *: the total flux density of the northern and 
central components. 



1999 



2004 



Comp. 


•S's.o 


•S's.o 


•58.4 


a 


® maj 




PA 




mJy 


mJy 


mJy 




mas 


mas 


deg 


N 




13 


9 


0.70 


1.07 


0.27 


28 


C 


830* 


644 


662 


-0.05 


0.50 


0.32 


50 


s 


364 


363 


269 


0.60 


0.70 


0.38 


40 



for 13 mJy. Component S is located at about 16.2 pc and 
in p. a. -146° with respect to component C, and accounts 
for 363 and 269 mJy at 5 and 8.4 GHz, respectively. 
Component N is located at about 20 pc and in p. a. 13° 
with respect to component C, and accounts for 13 and 9 
mJy at 5 and 8.4 GHz respectively. Both components S 
and N have a steep spectrum with afio = 0.60 ± 0.15 and 
Qsio = O-^O ± 0.30 respectively. Furthermore, with a flux 
density ratio Ss/^n = 28 - 30 , depending on the frequency, 
component N and S are very asymmetric also in their 
observed luminosity. 



3.3 Proper motion 

Decisive evidence supporting the idea that CSOs are radio 
sources in an early evolutionary stage has been provided 
by the measurements of the proper motion of the hot spots 
in compact CSOs (see e.g. iPolatidis fc Conwavl |2003| '1. To 
date, the CSO J0650-I-6001 lacks a proper measurement of 
expansion velocity. On the other hand, an apparent source 



contraction has been reported by lAkuior et all (|l996l ). rep- 
resenting a very uncommon case among radio sources. 
To derive a possible source expansion, we decided to mea- 
sure the distance between components S and N from the 
VLBA+EfTelsberg image at 5 GHz. Then we compared the 
result with what reported in IXu et all (|l995l ') and we found 
that the separation between these two components has in- 
creased by 0.28±0.13 mas in 13.5 years (i.e. the time interval 
between the two observations). We derive an apparent sep- 
aration velocity of Ws,a = 0.39c ± 0.18c, that corresponds to 
a kinematic age of tkin = 360 ±170 years. The accuracy of 
this separation velocity is limited since only two-epoch ob- 
servations are available with adequate resolution and at the 
same frequency. We do not cons ider the VLBA observation s 
carried out in 1999 and 2006 l|Orienti fc Dallacasal liooj ') . 
since the resolution was not adequate to reliably constrain 
the position of the northern component. 
To verify the apparent contraction reported bv lAkuior et al.l 
( 199d ). we determined the separation between the cen- 
tral and the southern components from 1999 to 2006, and 



we com pared our resu lts with those r eported by IXu et al 
(199^, iBritzen et aD (l2007h . and lOrienti fc Dallacasa 
((2003) (see Table [4|. Then we modelled the data with a 
linear fit and we obtain that the separation between com- 
ponents C and S has decreased of about 0.29±0.02 mas in 
15.6 years (i.e. the time interval spanned by the observa- 
tions), corresponding to an apparent "contraction" velocity 
'^c.a = 0.37c ± 0.02c (Fig. [3|. This suggests that component 
C likely enshrouds both the core region and the jet base 
which cannot be resolved with the spatial resolution (~0.6 
mas) provided by the observations at the highest frequency. 
The observed contraction is indicative of a knot in the jet 
that is moving towards the southern lobe and it is oriented 
with a small angle with respect to the line of sight. For this 
reason the flux density of component C is likely dominated 
by this knot whose emission is enhanced by boosting effects. 
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objects identified with galaxies, where boosting effects are 
thought to play a mi nor role. Nonetheless , in so me of the 
sources presented in iPolatidis fc Conwavl l|2003l l also jets 
have been detected and the proper motion of knots have 
been measured (for example in the sources J1945+7055 and 
J2355+4950). The target source J0650+6001 is optically 
identified with a quasar, and its peculiar characteristics, 
such as the radio emission dominated by the core compo- 
nent, the high flux density asymmetry between the outer 
components and their misalignment, may be due to orienta- 
tion effects that may modify the intrinsic source properties. 
To verify this hypothesis, we consider a beaming model 
which assumes that the source is intrinsically symmetric 
and the different flux density displayed by components 
N and S is due to Doppler boosting. In this context, the 
flux density ratio is related to both the jet speed and its 
orientation Q with respect to the line of sight: 



Figure 4. The light curve of J0650-I-6001 at 5 GHz. Blue crosses 
indicate the source flux density from VLA data (Table 2), while 
red triangles and green squares refer to VLBI flux density of com- 
ponent C and S respectively. The error bars refer to 3% and 5% 
of the VLA and VLBI flux density respectively (see Section 2). 
In addition to the VLBI data presented in this paper, we include 
those from Stanghellini et al. (1999), Xu et al. (1995), and Orienti 
& Dallacasa (2008). 



3.4 Flux density variability 

The light curve at 5 GHz of J0650-I-6001 clearly indicates a 
steady increment of the total flux density measured by the 
VLA over a period of nearly 18 years (see Section 2). The 
presence of flux density variability in a CSO is somewhat 
surprising, since young radio sources have b een considere d 
the least variable class of radio sources lO'Deal Il998l ). 
However, evidence of an increment in the flux density in 
the optically-thick part of the spectrum may be easily 
explained in terms of source expansion, assuming that the 
spectral peak is due to SS A, as in the case of J0650-I-6001 
l|Orienti fc Dallacasal [iooj ) . the source discussed here. As 
the source grows, the peak shifts towards lower frequencies, 
while the flux den sity in the optically-thick regime increases 
l|Pacholczvklll97d ). 

In Fig. |4] together with the total VLA flux density, we 
report also the VLB A flux density of components C and S. 
In the light curve of both the total source and component 
C, in addition to the steady increment of the flux density, 
it is also present a further contribution detectable in the 
1998 (VLA) and 1999 (VLBA) data. On the other hand, 
the light curve of component S is almost constant, with 
a possibly slight steady decrease, but without evidence of 
significant variability. This behaviour suggests that the 
total flux density is dominated by component C, that likely 
had an outburst close to the 1998-1999 observing epoch. 



4 DISCUSSION 

The determination of the properties of the population of 
young radio sources is mainly based on the analysis of 



Ss_ ^ / l-t-/3cosg \ 
Sn V 1 - /?cos6» ) 



3 + a 



(1) 



where we have assumed that component S is the one on the 
approaching side and P = v/c, where c is the speed of light. 
From VLBI data we found that R = 28 - 30, depending on 
the frequency considered, and a = 0.7 ± 0.3. With these 
values in Eq. [1] we find the possible combination of /3 and 
9 that is reported in Fig. [5] between dotted lines. We must 
note that although Eq. [1] usually applies to jet components, 
it can be used also to our purpose, since components N and 
S are hotspot dominated and an additional contribution 
from extended feature s like lobes is marginal sin ce their 
high magnetic field l|Orienti fc Dallacasal |2008| ) causes 
strong radiative losses at least at the frequencies considered 
here. 

Another way to derive the (/3-S) combination is from the 
measured separation velocity /3s, a between components N 
and S, by means of: 



2/3s,iSin6l 
r~Sf^cos2 



(2) 



Considering u^.a = 0.39c ± 0.18c as determined in Section 
3.3, the resulting possible {13-6) combinations are reported 
in Fig. [5] between solid lines. Comparing the areas obtained 
with the two different methods, we obtain an intrinsic 
separation velocity Ds,i = 0.43c ± 0.04c and a source 
orientation with respect to the line of sight of 12° and 28°. 
The comparison of the flux density of the outer components 
may be not always appropriate. In fact, given the relativistic 
time dilation, we are watching the components at different 
stages of their evolution. This is particularly important in 
the case of blazars whose jets are oriented with very small 
angles 6 < 5° to the line of sight, causing a quite rapid 
evolution due to boosting effects. On the other hand, this 
effect should not be relevant for J0650-I-6001 if we consider 
that the lobes/hotspots of this source are oriented at a 
relatively larger angle, as also supported by the lack of flux 
density variability in the southern component. 
Another way to determine the (/3-0) combination is based 
on the comparison of the distance of the two lobes with 
respect to the core. However, in this source the central 
component is likely dominated by a knot in the jet, thus 
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precluding the determination of the core position. The 
apparent contraction between component C and S is clearly 
indicating the presence of a compact component moving 
towards South with apparent speed iic,a = 0.37c ± 0.02c. 
If in Eq. [2] we consider the apparent separation velocity 
computed in Section 3.3 and the same orientation as derived 
from Fig. (5] we find that this knot has a mildly relativistic 
intrinsic velocity «c,i = 0.66c ± 0.03c. 

With the derived values we can compute the Doppler factor 
of the jet knot by means of: 



rc,i(l-/3c,iCose) ^ ' 

where Fc,! = — /S^j) is the Lorentz factor. 

If in Eq. Owe consider that the knot in the jet in component 
C is moving with Wc,i = 0.66c ± 0.03c, and it is oriented 
with an angle between 12° and 28° to the line of sight, we 
derive a Doppler factor D ~ 2, which is consistent with the 
peculiarly high observed flux density shown by the central 
component. In the case the whole flux density of component 
C is attributed to a knot in the jet, we can compute the 
intrinsic flux density S'lntr: 



ca 




Figure 5. The (e,/3)-plane for J0650+6001, as derived from 
the northern and southern components. The solid lines repre- 
sent the limits derived from the apparent separation velocity 
(/3s, a = 0.39 ±0.18), while the dotted lines refer to limits inferred 
from the flux density ratio. 



Sintr = Sobs^-''+"' (4) 

where Sobs is the observed flux density. If in Eq. 2] we 
consider the observed flux density reported in Table O 
D ~ 2 as derived above, and a — 0.7, we flnd that the 
intrinsic flux density of component C, i.e. core+jet region, 
is about 50 mjy at both 5 and 8.4 GHz, in good agreement 
with what expected for the central component of radio 
sources. 

Another explanation for the misaligned source structure 
and its high flux density asymmetry may be caused by 
an inhomogeneous external environment. In fact, in the 
case the two jets/hotspots experience different ambient 
properties on the two sides of the radio source, the side 
hitting the most dense medium has its expansion slowed 
down and its radio emission enhanced, since the energy 
losses due to adiabatic expansion are smaller with respect to 
the other side. This ends up with the brighter and smaller 
radio component lying closer to the core. If this is the case, 
the model developed assuming beaming effects becomes 
meaningless, and the source may thus be oriented even at 
larger angles to the line of sight, making the estimation of 
the intrinsic velocity an upper limit. 



a source kinematic age of 360±170 years. On the other 
hand, the distance between the central component and 
the southern lobe has decreased with a mean apparent 
contraction speed «c,a = 0.37c ± 0.02c. To explain these 
peculiar characteristics, we discussed an interpretation 
based on relativistic beaming effects. From this analysis we 
drew a picture in which the source is oriented with an angle 
between 12° and 28° to the line of sight with the southern 
component being the approaching side of the source. The 
intrinsic separation velocity is Vsa — 0.43c ± 0.04c. The 
apparent contraction observed between the central and 
the southern components may be explained in terms of a 
knot in the jet that dominates the radio emission of the 
central component, that is moving towards the southern 
component with a mildly relativistic intrinsic velocity 
Dc,i = 0.66c ± 0.03c, and with a Doppler factor D = 2, 
which causes a moderate boost of its flux density. The 
monitoring of the flux density variability and the study of 
the component separations on a longer time scale will be 
crucial to better determine the nature of this source. 



5 CONCLUSIONS 

We have presented multi-epoch VLBI and VLA observa- 
tions of the radio source J0650+6001. The source shows 
a peculiar triple structure, where the outer components 
are misaligned with respect to the central flat-spectrum 
region, and show a high flux density asymmetry. The 
comparison between multi-epoch high spatial resolution 
observations indicated that the farthest components are 
increasing their distance with an apparent separation 
velocity Vs,a = 0.39c ± 0.18c. Assuming that this velocity 
is representative of the mean separation speed, we obtain 
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